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The angular distributions and the partial branching fraction of the decay — > K*'^ fi'^ fi'^ are 
studied using an integrated luminosity of 0.37 fb~^ of data collected with the LHCb detector. The 
forward-backward asymmetry of the muons, Afb, the fraction of longitudinal polarisation, Fl, and 
the partial branching fraction, dB/dq^, are determined as a function of the dimuon invariant mass. 
The measurements are in good agreement with the Standard Model predictions and are the most 
precise to date. In the dimuon invariant mass squared range 1.00 - e.OOGcVVc", the results are 
^FB = -0.06 l^ i|± 0.04, Fl = 0.55±0.10±0.03 and dB/dq^ = (0.42 ± 0.06 ± 0.03) x lO'^c''/ GeV^ 
In each case, the first error is statistical and the second systematic. 
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The process — > is a flavour changing neu- 

tral current decay. In the Standard Model (SM) such 
decays are suppressed, as they can only proceed via loop 
processes involving electrowcak penguin or box diagrams. 
As-yet undiscovered particles could give additional con- 
tributions with comparable amplitudes, and the decay is 
therefore a sensitive probe of new phenomena. A num- 
ber of angular observables in — > decays 
can be theoretically predicted with good control of the 
relevant form factor uncertainties. These include the 
forward-backward asymmetry of the muons, ApB, and 
the fraction of longitudinal polarisation, Fl, as functions 
of the dimuon invariant mass squared, [1]. These ob- 
servables have previously been measured by the BaBar, 
Belle, and CDF experiments [2]. A more precise de- 
termination of A-p-Q is of particular interest as, in the 
1.00 < < 6.00GeV^/c'* region, previous measurements 
favour an asymmetry with the opposite sign to that ex- 
pected in the SM. If confirmed, this would be an unequiv- 
ocal sign of phenomena not described by the SM. This 
letter presents the most precise measurements of Ape, 
and the partial branching fraction, dS/dg^, to date. 
The data used for this analysis were taken with the LHCb 
detector at CERN during 2011 and correspond to an inte- 
grated luminosity of 0.37 fb~^. The is reconstructed 
through its decay into the K~^Tr~ final state. 

The LHCb detector [3] is a single-arm spectrometer 
designed to study 6-hadron decays. A silicon strip ver- 
tex detector positioned around the interaction region is 
used to measure the trajectory of charged particles and 
allows the reconstruction of the primary proton-proton 
interactions and the displaced secondary vertices charac- 
teristic of B-meson decays. A dipole magnetic field and 
further charged particle tracking stations allow momenta 
in the range 5 < p < 100 GeV/c to be determined with 
a precision of 6p/p = 0.4-0.6%. The experiment has an 
acceptance for charged particles with pseudorapidity be- 
tween 2 and 5. Two ring imaging Cherenkov (RICH) de- 
tectors allow kaons to be separated from pions or muons 
over a momentum range 2 < p < 100 GeV/ c. Muons are 
identified on the basis of the number of hits in detectors 
interleaved with an iron muon filter. 

The i?° —7- K*^ijl^ii~ angular distribution is governed 
by six (j^-dcpendcnt transversity amplitudes. The decay 
can be described by q^ and the three angles 0;, Ok-, 
For the (S"), Si is the angle between the /i+ and 
the opposite of the {B^) direction in the dimuon rest 
frame, 9k the angle between the kaon and the direction 
opposite to the B meson in the K*^ rest frame, and </> the 
angle between the jJL'^ and if+Tr" decay planes in the 
B rest frame. The inclusion of charge conjugate modes is 
implied throughout this letter. At a given g^, neglecting 
the muon mass, the normalised partial differential width 
integrated over 6k and (j) is 



^'5^ , =^FL(l-COS^gO + 

T d COS 01 Aq^ 4 

1(1 - i^L)(l + COS^ Oi) + ^FB COS^, 
o 

and integrated over 6i and <f) it is 



(1) 



1 



d^F 



9k + 



3 2 

Fdcos0xdg2 - 2^"^°^ ' 



-(1-Fl)(1-cos2 



(2) 



of 5°^ K*°u+ij,- decays to the yield of B°^ J/ipK*° 
11+ H~. 



These expressions do not include any broad S-wave con- 
tribution to the B^ K+TT~ fi^ decay and any con- 
tribution from low mass tails of higher K*'^ resonances. 
These contributions are assumed to be small and are ne- 
glected in the rest of the analysis. 

Signal candidates are isolated from the background us- 
ing a set of selection criteria which are detailed below. An 
event-by-event weight is then used to correct for the bias 
induced by the reconstruction, trigger and selection crite- 
ria. In order to extract ApB and Fl, simultaneous fits are 
made to the K'^'k~ ^i^ ^~ invariant mass distribution and 
the angular distributions. The partial branching fraction 
is measured by comparing the efficiency corrected yield 

where J/ip —. 

Candidate K*^n'+ ^~ events are first required to 

pass a hardware trigger which selects muons with a trans- 
verse momentum, > 1.48GcV/c. In the subsequent 
software trigger, at least one of the final state particles is 
required to have both px > 0.8 GeV/c and impact param- 
eter > 100 fiva with respect to all of the primary proton- 
proton interaction vertices in the event [4]. Finally, the 
tracks of two or more of the final state particles arc re- 
quired to form a vertex which is significantly displaced 
from the primary vertices in the event [5]. 

In the final event selection, candidates 
with i4r+7r^/i+/i^ invariant mass in the range 
5100 < m^+TT-Ai+M- < 5600 MeV/c^ and A'+tt" in- 
variant mass in the range 792 < < 992MeV/c^ 
are accepted. Two types of backgrounds arc then 
considered: combinatorial backgrounds, whcirc the 
particles selected do not come from a single 6-hadron 
decay; and peaking backgrounds, where a single de- 
cay is selected but with some of the particle types 
mis-identified. In addition, the decays B^ ^ J/ijjK*^ 
and S"^ Vj(25)A:*", where J/ip ,^{2S) ^ n+ n- , are 
removed by rejecting events with dimuon invariant mass, 
m^+^-, in the range 2946 < to^+^- < 3176MeV/c^ or 
3586 < m^+f,- < 3776MeV/c2. 

The combinatorial background, which is smoothly 
distributed in the reconstructed Tr~ iJ,~^ fi~ invariant 
mass, is reduced using a Boosted Decision Tree (BDT). 
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The BDT uses information about the event kinematics, 
vertex and track quality, impact parameter and particle 
identification information from the RICH and muon de- 
tectors. The variables that are used in the BDT are cho- 
sen so as to induce the minimum possible distortion in the 
angular and distributions. For example, no additional 
requirement is made on the of both of the muons as, 
at low , this would remove a large proportion of events 
with I cos6'/| ~ 1. The BDT is trained entirely on data, 
using samples that are independent of that which is used 
to make the measurements: triggered and fully recon- 
structed i?" — > J/ipK*^ events are used as a proxy for the 
signal decay, and events from the upper — ^ K*^ fX^ 
mass sideband (5350 < mj^+j^- ^+^- < 5600MeV/c^) are 
used as a background sample. The lower mass sideband 
is not used, as it contains background events formed from 
partially reconstructed B decays. These events make a 
negligible contribution in the signal region and have prop- 
erties different from the combinatorial background which 
is the dominant background in this region. 

A cut is made on the BDT output in order to optimise 
the sensitivity to ^fb averaged over all q^. The selected 
sample has a signal-to-background ratio of three to one. 

Peaking backgrounds from — > </)/i^"/i^ (where 
(j) K+K-), B° J/iPK*° and iIj{2S)K*° are 

considered and reduced with a set of vetoes. In each 
case, for the decay to be a potential signal candidate, at 
least one particle needs to be misidentified. For example, 
B"— > J/ip K*'^ events where a kaon or pion is swapped for 
one of the muons, peak around the nominal B'^ mass and 
evade the J/i/j veto described above. Vetoes for each of 
these backgrounds are formed by changing the relevant 
particle mass hypotheses and recomputing the invariant 
masses, and by making use of the particle identification 
information. In order to avoid having a strongly peak- 
ing contribution to the cos 9k angular distribution in the 
upper mass sideband, B^ — > fi^ ^~ candidates are 
removed. Events with A'+/i+/^^ invariant mass within 
60 MeV/c^ of the nominal B'^ mass are rejected. The ve- 
toes for all of these peaking backgrounds remove a neg- 
ligible amount of signal. 

After the application of the BDT cut and the above ve- 
toes, a fit is made to the K^tt~ fi'^ fj,^ invariant mass dis- 
tribution in the entire accepted mass range (see Fig. 1). A 
double-Gaussian distribution is used for the signal mass 
shape and an exponential function for the background. 
The signal shape is fixed from data using a fit to the 
— > J/t/jK*'^ mass peak. In the full q^ range, in a 
signal mass window of ±50MeV/c^ (±2.5cr) around the 
measured i?° mass, the fit gives an estimate of 337 ±21 
signal events with a background of 97 ± 6 events. 

The residual peaking background is estimated using 
simulated events. As detailed below, the accuracy of the 
simulation is verified by comparing the particle (mis-) 
identification probabilities with those derived from con- 
trol channels selected from the data. The residual peak- 




FIG. 1. K^TT^ jj.'^ fi'^ invariant mass distribution after the 
application of the full selection as data points with the fit 
overlaid. The signal component is the green (light) line, the 
background the red (dashed) line and the full distribution the 
blue (dark) line. 

ing backgrounds are reduced to a level of 6.1 events, i.e. 
1.8% of the 337 observed signal events. The backgrounds 
from 5° — > and B'^ — > J/iJjK*° decays do not 

give rise to any forward-backward asymmetry and are 
ignored. However, in addition to the above backgrounds, 
— >■ decays with the kaon and pion swapped 

give rise to a 0.7% contribution. The change in the sign 
of the particle which is taken to be the kaon results in 
a B^ {B'^) being reconstructed as a B^ {B°), therefore 
changing the sign of A^b for the candidate. This misiden- 
tification is accounted for in the fit for the angular ob- 
servables. 

The selected 5°— >■ K*'^ jjL^ candidates are weighted 
in order to correct for the effects of the reconstruction, 
trigger and selection. The weights are derived from sim- 
ulated B^ — )■ K*^fi~^fj,~ events and are normalised such 
that the average weight is one. In order to be indepen- 
dent of the physics model used in the simulation, the 
weights are computed based on cos^i^-, cos^; and q^ on 
an event-by-event basis. The variation of detector effi- 
ciency with the (j) angle is small and ignoring this varia- 
tion does not bias the measurements. Only events with 
0.10 < g2 < 19.00 GeVVc'^ are analysed. 

Owing to the relatively unbiased selection, 89% of 
events have weights between 0.7 and 1.3, and only 3% 
of events have a weight above 2. The distortions in the 
distributions of cos 6k, cos 9i and q^ that are induced 
originate from two main sources. Firstly, in order to pass 
through the iron muon filter and give hits in the muon 
stations, tracks must have at least 3GeV/c momentum. 
At low q"^ this removes events with | cos 9i\ ^ 1. This ef- 
fect stems from the geometry of the LHCb detector and is 
therefore relatively easy to model. Secondly, events with 
cos 9k ~ 1, and hence a slow pion, are removed both 
by the pion reconstruction and by the impact parameter 
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requirements used in the trigger and BDT selection. 

A number of control samples are used to verify the 
simulation quality and to correct for differences with 
respect to the data. The reproduction of the i?° mo- 
mentum and pseudorapidity distributions is verified us- 
ing i?" — !• J/ijjK*^ decays. These decays are also used 
to check that the simulation reproduces the measured 
properties of selected events. The hadron and muon 
(mis-)identification probabilities are adjusted using de- 
cays where the tested particle type can be determined 
without the use of the particle identification algorithms. 
A tag and probe approach with J/tp — > fJ.'^ fJ.~ decays is 
used to isolate a clean sample of genuine muons. The 
decay D*+ D°7r+, where 1?° K~tt^, is used to 
give an unambiguous source of kaons and pious. The 
statistical precision with which it is possible to make the 
data/simulation comparison gives rise to a systematic un- 
certainty in the weights which is evaluated below. 

The observables A-pB and Fl are extracted in bins of 

. In each bin, a simultaneous fit to the X+Tr^/i+zi^ 
invariant mass distribution and the cos Ok and cos 9i dis- 
tributions is performed. The angular distributions are 
fitted in both the signal mass window and in the upper 
mass sideband which determines the background param- 
eters. The angular distributions for the signal are given 
by Eqs. 1 and 2 and a second order polynomial in cos 9k 
and in cos 9i is used for the background. 

In order to obtain a positive probability density func- 
tion over the entire angular range, Eqs. 1 and 2 imply 
that the conditions |Afb| < |(1 - -Fl) and < Fl < 1 
must be satisfied. To account for this, the maximum like- 
lihood values for ^fb and Fl are extracted by performing 
a profile-likelihood scan over the allowed range. The un- 
certainty on the central value of ApB and Fl is calculated 
by integrating the probability density extracted from the 
likelihood, assuming a flat prior in ApB and Fl, inside 
the allowed range. This gives an (asymmetric) 68% con- 
fidence interval. 

The partial branching fraction is measured in each 
of the bins from a fit to the efficiency corrected 
J^+TT^/i+Zi^ mass spectrum. The efficiencies are deter- 
mined relative to the J/ipK*^ decay which is used 
as a normalisation mode. 

The event weighting and fitting procedure is vali- 
dated by fitting the angular distribution of — ^ J/tpK*^ 
events, where the physics parameters are known from 
previous measurements [6]. The product of the 
B'^ ~> J/tpK*'-' and J/tjj — > /i^/i^ branching fractions 
is ^ 75 times larger than the branching fraction of 
S*^— > K*^ IJ& , allowing a precise test of the procedure 
to be made. Fitting the — > J/ip FT*" angular distribu- 
tion, weighted according to the event-by-event procedure 
described above, yields values for Fl and A^b in good 
agreement with those found previously. 

For -J> K*°fj,+n-, the fit results for Ape, and 
dB/dq^ are shown in Fig. 2 and are tabulated together 
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FIG. 2. Afb, Fl and dB/dq^ as a function of . The SM pre- 
diction is given by the cyan (light) band, and this prediction 
rate-averaged across the q^ bins is indicated by the purple 
(dark) regions. No SM prediction is shown for the region be- 
tween the two regimes in which the theoretical calculations 
are made (see text). 

with the signal and background yields in Table I. The fit 
projections are given in the appendix. Signal candidates 
are observed in each q^ bin with more than 5a signif- 
icance. The compatibility of the fits and the data are 
assessed using a binned test and all fits are found to 
be of good quality. The measurements in all three quan- 
tities are more precise than those of previous experiments 
and are in good agreement with the SM predictions. The 
predictions are taken from Ref. [7]. In the low q^ region 
they rely on the factorisation approach [8] , which loses ac- 
curacy when approaching the J/^p resonance; in the high 
q^ region, an operator product expansion in the inverse 
6-quark mass, l/rrib, and in l/\/<f is used [9], which is 
only valid above the open charm threshold. In both re- 
gions the form factor calculations are taken from Ref. [10] 
and a dimensional estimate is made on the uncertainty 
from expansion corrections [11]. 

In the 1.00 < < G.OOGeVVc'^ region, the fit gives 
Afb = -Omto ll ± 0.04, Fl = 0.55 ± 0.10 ± 0.03 and 
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dB/dq^ = (0.42 ±0.06 ±0.03) x 10-V/GeV^ where the 
first error is statistical and the second systematic. The 
theoretical predictions in the same range are ^fb = 
-0.04±0.03, Fl = 0.74;°;°? and dS/dg^ = (0.50 x 
lO-'^c^/GeV^ The LHCb Afb measurement is a factor 
1.5 — 2.0 more precise than previous measurements from 
the Belle, CDF and BaBar collaborations [2] which are, 
respectively, ApB = 0.261°;^^ ± 0.07, Apg = 0.29^°;^° ± 
0.07 and, for < 6.25 GeVVc'^, ^fb = 0.24^°:^^ ± 0.05. 
The positive value of ApB preferred in the 1.00 < q^ < 
6.00 GeV^/c'^ range in these previous measurements is not 
favoured by the LHCb data. The previous measurements 
of i^L in the same regions are i^L = 0.67 ± 0.23 ± 0.05 
(Belle), Fl = OMta ^^i ± 0.08 (CDF) and Fl = 0.35 ± 
0.16 ± 0.04 (BaBar). These are in good agreement with 
the LHCb resuh. 



For the determination of Afb and Fl, the dominant 
systematic uncertainties arise from the event-by-event 
weights which are extracted from simulated events, and 
from the model used to describe the angular distribution 
of the background. The uncertainty on the event-by- 
event weights is evaluated by fluctuating these weights 
within their statistical uncertainties and repeating the 
fitting procedure. The uncertainty from the background 
model which is used is estimated by changing this model 
to one which uses binned templates from the upper mass 
sideband rather than a polynomial parameterisation. 

The dominant systematic errors for the determination 
of dB/dq^ arise from the uncertainties on the particle 
identification and track reconstruction efficiencies. These 
efficiencies are extracted from control channels and are 
limited by the relevant sample sizes. The systematic 
uncertainty is estimated by fluctuating the efficiencies 
within the relevant uncertainties and repeating the fit- 
ting procedure. An additional systematic uncertainty of 
~ 4% arises from the uncertainty in the i?" — > J/tpK*'^ 
and J/ip — >■ branching fractions [12]. 

The total systematic error on each of ^fb and Fl 
(dB/dq'^) is typically - 30% (50%) of the statistical error, 
and hence adds ~ 4% (~ 11%) to the total uncertainty. 

In summary, using 0.37 fb~ of data taken with the 
LHCb detector during 2011, Ape, Fl and dB/dq^ have 
been determined for the decay 5°— K*^fi^^~. These 
are the most precise measurements of these quantities 
to-date. All three observables show good agreement with 
the SM predictions. 
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TABLE I. Central values with statistical and systematic uncertainties for Afb, and dB/dq as a function of q^. The 
B° — > K*'^ jX^ signal and background yields in the ±50MeV/c^ signal mass window with their statistical uncertainties are 
also indicated, together with the statistical significance of the signal peak that is observed. 









^FB 




Fl 


dB/dq^ 


Signal 


Background 


Significance 


(GeVVc") 










(xlO-'' cVGeV^) 


yield 


yield 


{<y) 


0.10 < q^ 


< 2.00 


-0.15 ±0.20 ±0.06 


0.00 ±0-02 


0.61 ±0.12 ±0.06 


48.6 ±8.1 


16.2 ±2.3 


8.6 


2.00 < q"^ 


< 4.30 


0.05 


+ 0.16 
- 0.20 


±0.04 


0.77± 0.15 ±0.03 


0.34 ±0.09 ±0.02 


26.5 ±6.5 


15.7± 2.2 


5.4 


4.30 < q^ 


< 8.68 


0.27 


+ 0.06 
- 0.08 


±0.02 


0.60 in? ±0.01 


0.69 ±0.08 ±0.05 


104.7 ± 11.9 


31.7±3.3 


12.4 


10.09 < 


< 12.86 


0.27 


+ 0.11 
- 0.13 


±0.02 


0.41 ±0.11 ±0.03 


0.55 ± 0.09 ± 0.07 


62.2 ± 9.2 


20.4 ± 2.6 


9.6 


14.18 < 


< 16.00 


0.47 


+ 0.06 
- 0.08 


±0.03 


0.37 ±0.09 ±0.05 


0.63 ±0.11 ±0.05 


44.2 ± 7.0 


4.2 ± 1.3 


10.2 


16.00 < q^ 


< 19.00 


0.16 


+ 0.11 
-0.13 


±0.06 


0.26 lg;o8 ±0.03 


0.50 ±0.08 ±0.05 


53.4 ±8.1 


7.0 ± 1.7 


9.8 


1.00 < q^ 


< 6.00 


-0.06 


+ 0.13 
- 0.14 


±0.04 


0.55 ±0.10 ±0.03 


0.42 ± 0.06 ± 0.03 


76.5 ± 10.6 


33.1 ±3.2 


9.9 



Appendix 

The following fit projections are published as EPAPS material. 
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FIG. 3. Fit projections for rriKTv^fj,, cos9i and cos6'x for the bins: 0.10 < < 2.00, 2.00 < < 4.30 and 4.30 < q^ < 
8.68GeVVc^ 
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FIG. 4. Fit projections for m^TrMM, cosdi and cos Ok for the bins: 10.09 < < 12.86, 14.18 < < 16.00 and 16.00 < < 
19.00 GeVVc*. 



